clear all;
it = 22;
act_name = ['ACTMAT',num2str(it),'.rosmat'];
last_name = ['LASTMAT',num2str(it),'.rosmat'];
obj_name = ['OBJ',num2str(it),'.rosmat'];
dref_name = ['DREF',num2str(it),'.rosmat'];
vref_name = ['aVREF',num2str(it),'.rosmat'];
Mact = load(act_name);
Mlast = load(last_name);
Mobj = load(obj_name);
Mdref = load(dref_name);
Mvref = load(vref_name);
% subplot(141)
% pcolor(Mact);
% shading interp
% subplot(142)
% pcolor(Mlast)
% shading interp
figure(1)
subplot(121)
pcolor(Mobj);
shading flat
max_obj = max(max(Mobj));
title(num2str(max_obj));
subplot(122)
pcolor(abs(Mact-Mlast));
shading flat
figure(2)
Mdif = Mact-Mlast;
Mdif(find(Mdif<0.0125))=0;
subplot(121)
spy(Mlast-Mact)
subplot(122)
spy(Mdif)
figure(3)
subplot(141)
pcolor(Mact);
shading interp
subplot(142)
pcolor(Mlast)
shading interp
subplot(143)
pcolor(Mdref)
shading interp
subplot(144)
pcolor(Mvref)
shading interp
